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Low frequency acoustic wave generation is studied taking into account the effect of stratification, 
inhomogeneity of background velocity profile and temperature fluctuations, ft is shown that for the 
typical parameters of convective storms the dipole radiation related to temperature inhomogeneities 
is at least of the same order as radiation of Lighthill's quadrupole source. It is also shown that 
the source related to stratification could have valuable contribution whereas some other sources are 
shown to be inefficient. 
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I. INTRODUCTION 

It is long known that strong convective storms, such 
as supercell thunderstorms are powerful sources of in- 
frasound Q, 1^ 1^. Detailed observations of convective 
storm generated infrasound 0, Q provides that at least 
two different group of infrasonic signals could be identi- 
fied. The first group, with characteristic period about 1 
s, has been found to be in strong connection with proto- 
tornadic structures, funnel clouds and tornadoes. Based 
on coincident radar measurements of tornados, which 
show strong relationship between funnel diameter and 
infrasound frequency, it is usually supposed that these in- 
frasound waves are generated by radial vibrations of the 
funnel core Q. The second group of infrasound signals 
has the periods from 2 to 60 s. Usually the emission ap- 
pears about 1 hour before observation of tornado. These 
waves are not related with tornado itself and are caused 
by convective processes that precedes tornado formation. 
The acoustic power radiated by convective storm system 
could be as high as 10'' watts (7j. Although several rea- 
sonable mechanisms have been suggested to explain this 
acoustic radiation, thephysical mechanism of the process 
remains unexplained 0, Ia.l3- 

Broad and smooth spectrum of the observed infra- 
sound radiation indicates that turbulence is one of the 
promising sources of the radiation. Lighthill's acoustic 
analogy Q represents the basis for understanding of the 
generation of sound by turbulent flows. In this approach 
the flow is assumed to be known and the sound field is 
calculated as a small by-product of the flow. According 
to this theory in the case of uniform background ther- 
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modynamic parameters interaction of turbulent vortices 
provides quadrupole source of sound. The acoustic power 
of the source was estimated by Proudman 1^. But us- 
age of this estimation for the infrasound radiation from 
convective storms usually leads to the underestimation 
of the acoustic power 0, Q ■ 

It is also well known that any kind of inhomogeneity of 
the background flow, such as stratification, shear of the 
background velocity profile and temperature fluctuations 
leads to appearance of additional (mainly dipole) sources 
of sound jlO,, JL IJ . Resent developments in this direction 
led to the formulation of the generalized acoustic analogy 
that implies: (i) dividing the the flow variables into 
their mean and fluctuating parts; (ii) substracting out 
the equation for the mean flow; (iii) collecting all the 
linear terms on one side of equations and the nonlinear 
terms on the other side; (iv) treating the latter terms as 
as known source of sound. 

In the presented paper we consider acoustic radiation 
from turbulent convection taking into account mentioned 
above effects of stratification, inhomogeneity of velocity 
profile and temperature fiuctuations. Performed anal- 
ysis shows that for the typical parameters of supercell 
storms the dipole radiation related to temperature inho- 
mogeneities is at least of the same order as radiation of 
Lighthill's quadrupole source. It is also shown that the 
source related to stratification could have valuable con- 
tribution whereas some other sources are shown to be 
very inefficient. 

The paper is organized as follows: simple background 
flow model is constructed in Sec. ^] Linear equation gov- 
erning the propagation of sound in the flow is obtained 
in Sec. IIIII Sources of acoustic radiation are obtained 
and analyzed in Sec. IIVI Summary is given is Sec. 
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II. 



BACKGROUND FLOW MODEL 



In this section we construct simplified model for back- 
ground updraft flow. Suppose in the stratified dry atmo- 
sphere there exists some region (see Fig. 1) with supply 
of relatively hot air at some background level and a sink 
at some height h. Assume pressure, density and temper- 
ature fields outside this region are Pq{z), Pq{z) and Tq{z) 
respectively. In the case of isentropic atmosphere Q 
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for isentropic atmosphere and 



pi(z) = po.^exp(--^) 



(10) 



(11) 



(12) 



where H = RToh/g is stratification length scale; Pofc 
and Tof, are pressure and temperature at the background 
level; z is vertical coordinate; g is gravity acceleration 
and 7 = Cp/cv is adiabatic index. 

In the case of isothermal atmosphere instead of Eq. 
we have 
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Suppose the temperature of the supplied hot air is 
Tib > Tob- The continuity, Euler, heat and state equa- 
tions governing the stationary adiabatic updraft motion 
in the region are: 
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for isothermal atmosphere. 
Eqs. © and yields 
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where L is horizontal length scale of region 1 . Therefore 
if iJ ^ L we conclude that 
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where Vi = V + V^; V(z) = [0,0, ^(z)] is up- 
draft motion velocity; z is vertical coordinate; V£)(r) = 
[a{z)x, a{z)y, 0] is horizontal velocity associated with di- 
vergence of updraft motion; and pi(z), Pi(z), Ti(z) are 
density, pressure and temperature respectively. 

The boundary condition at the boundary of regions 1 
and 2 implies 



Pi(z) = Po(z). 
Eq. (O provides for a{z) 
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For other variables straightforward calculations yield 
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Considered model of background flow do not include 
horizontal rotation of the updraft flow as well as horizon- 
tal wind with vertical shear. Existence of the latest one is 
known to be necessary for the formation of the strongest 
convective storm systems, supercell storms 0,0. How- 
ever, as it will be shown below none of these flows can 
have valuable direct influence on the infrasound genera- 
tion. 



III. LINEAR OPERATOR FOR ACOUSTIC 
WAVES 

Next step of the study in the framework of mentioned 
above generalized acoustic analogy is derivation of the 
linear equation for acoustic wave propagation. Taking 
into account Eq. H16|l linearized continuity, Euler and 
state equations can be written as 
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FIG. 1: Schematic illustration of the considered updraft flow 
model. 



perturbations are very inefficient sources of sound. The 
acoustic power is proportional to e~^/^*^ and q—^^/^'^ 
for instability waves and continuous spectrum perturba- 
tions respectively ^3? ]. In the last expression S is the 
ratio of length scales of energy containing vortices and 
background velocity inhomogeneity {V/dzV). In the case 
of supercell thunderstorm M ~ 0.1 — 0.15 and 6 ^ 10~^, 
therefore both linear mechanisms have negligible acoustic 
output and attention should be payed to sources of sound 
related to nonlinear terms and entropy fluctuations that 
will be studied in the next section. 

Propagation of infrasound in the atmosphere was in- 
tensively studied by different authors ( see, e.g., [Tsl 
and references therein) and will not be considered in this 
paper. 

IV. SOURCES OF INFRASOUND 
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where v, p and p are perturbations of velocity, density 
and pressure respectively; Dt = dt + VdzV is convective 
derivative; and Dt = Dt + a. 

Combining Eqs. H17|) - (|21|) for perturbation of total en- 
thalpy 



Pi 

we obtain following equation 
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For obtaining the sources of acoustic waves generated 
by turbulent convection we follow standard technique of 
acoustic analogy. Dividing the flow variables into their 
mean and fluctuating parts keeping the nonlinear terms 
in Eqs. H17|l - H20|) . and noting that in general in the flow 
there can exist entropy fluctuations, so replacing Eq. 
by 
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where s is entropy fluctuation, after straightforward cal- 
culations we obtain following equation 



LQ^Si + Ss + Sv + St, 



where 

Si = dt[dxdi{viVx)+dydi{viVy)]+Dt 



dz - XjDt 



(26) 

diiviVz), 
(27) 



where L is so-called third order linear wave operator 



L = DtDt\Dt 
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dzPidz 
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and Cs = (7Pi/pi)^/^ is sound velocity and V'^ = d^+dy. 

Operators similar to L are well known in jet noise the- 
ory [O, lial . Obtained linear equation describes not only 
the acoustic wave propagation but it also describes the 
instability wave solutions that are usually associated with 
large scale turbulent structures and continuous spectrum 
solutions that are related to " fine-grained" turbulent mo- 
tions |l2l IT5I I . In the presence of any kind of inhomogene- 
ity (in the case under consideration all the background 
variables depend on vertical coordinate z) linear coupling 
between these perturbations is possible, and in principle 
acoustic waves can be generated by both instability waves 
and continuous spectrum perturbations. But in the case 
of low Mach number (M = V/cg 1) fiows both kind of 
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are source terms. In the derivation of Eq. (|26|l all the 
nonlinear terms containing pvi has been omitted, due to 
the fact that these terms are known to be much weaker 
sources of sound then the terms containing ViVj in the 
low Mach number fiows 16]. 

All the sources presented in Eq. H26|l are well known in 
aeroacoustics. Consequently, we do not repeat calcula- 
tions of acoustic power (i.e., acoustic energy generated in 
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a unit time) of each source term and adopt well known re- 
sults 0,0,0,113 for the problem under consideration. 
The main idea of mentioned methods of acoustic power 
calculations is the following: in the low Mach number 
flows convective propagation effects have minor influence 
on sound generation process and therefore all the terms 
containing V on the left hand side of Eq. (|2(j|l can be ne- 
glected. This circumstance allows one to obtain simple 
estimations of acoustic power of the sources. 

In the absence of background flow {V — 0) the flrst 
term Si reduces to the Lighthill's quadrupole source 0| . 
Its acoustic power 
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where v is rms of turbulent velocity fluctuations and F is 
total volume of turbulent flow. This expression was used 
in ref. Q to estimate acoustic power of supercell storm. 

The term Ss is dipole source caused by the stratifi- 
cation. The influence of stratification on acoustic wave 
generation has been studied intensively in the context of 
solar physics 0, |23| ■ Related acoustic power 



F = 



Ni, 



(32) 



where lo = v /I is the characteristic frequency of emitted 
acoustic waves. Substituting Cg ~ 340 m/s, H = 
m and u; ~ 0.1 s~^, we obtain Ng ^ O.lNi, i.e., the 
power of this source is weaker then acoustic power of 
Lighthill's quadrupole source, but in principle considered 
dipole source could have valuable contribution to the to- 
tal acoustic power. 

The source term Sv is related to inhoniogeneity of up- 
draft velocity. It is dipole source similar to Sg- Corre- 
sponding acoustic power 



(33) 



is negligibly small in comparison with Ng for low Mach 
number flows (for supercell storms M 0.1 — 0.15). 

The last term St is caused by entropy fluctuations. 
The most strong acoustic source of this kind is so-called 
thermo acoustical source 0| that is related to small scale 
inhomogeneities of background density (and therefore 
temperature) fluctuations, that produce dipole source. 
The physics of this kind of acoustic radiation is the fol- 
lowing: "hot spots" or " entropy inhomogeneities" behave 
as scattering centers at which dynamic pressure fluctu- 
ations are converted directly into sound. The acoustic 
power 
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where AT denotes the rms of temperature fluctuations. 

Taking for supercell storm v ~ (3—10) m/s; AT ~ 
(3 - 10)°K and T = 270° K we see that for the typical 
parameters of supercell storm the dipole radiation related 
to temperature inhomogeneities is at least of the same 
order as radiation of Lighthill's quadrupole source. 

In the considered model no attention was payed to the 
horizontal wind with vertical shear and rotation of the 
updraft flow. These assumptions allow us to obtain rel- 
atively simple equations. These flows could not make 
valuable changes in the acoustic output estimations. In- 
deed, consideration of the wind and rotation can lead 
to appearance in the right hand side of Eq. H2t)|l ad- 
ditional dipole terms similar to Sy that has negligible 
acoustic radiation, and also additional quadrupole source 
term known as " shear noise" . Acoustic power of the 
later one is also much weaker compared to Lighthill's 
quadrupole source in low Mach number flows. 

In the presented paper we study dry atmosphere and 
therefore influence of humidity on the sound generation 
in the convective storms has not been considered. We 
intend to study this influence in the future. 



V. SUMMARY 

The simple model of the background flow has been con- 
structed and the problem of acoustic radiation from tur- 
bulent convection taking into account the effects of strat- 
ification, inhomogeneity of velocity profile and tempera- 
ture fluctuations has been considered in the framework 
of acoustic analogy. Performed analysis shows that for 
the typical parameters of supercell storms the dipole ra- 
diation related to temperature inhomogeneities is at least 
of the same order as radiation of Lighthill's quadrupole 
source. It is also shown that the source related to strati- 
fication could have valuable contribution whereas dipole 
and quadrupole sources related to the inhomogeneity of 
background velocity profile are shown to be very ineffi- 
cient. 
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